Shallow Coulomb Gap and Weak Level-Level Correlations in a Deeply Insulating 

Electron System with Interactions 
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We consider a system of two-dimensional electrons strongly localized by disorder. Interactions 
create a gap in the average tunneling density of states v{E) at energies, E, close to the Fermi level. 
We derive a system of self-consistent equations for the correlators of local densities of occupied and 
empty states. When either the interactions are screened by the gate or the temperature, T, is high 
enough, so that the Coulomb gap is shallow, the perturbative solution of the system yields analytical 
expressions for the interaction-induced correction 8u(E) and the level-level correlator. We show that 
even with short-range interactions, 5v(E) exhibits a singular energy dependence at T = 0. We also 
demonstrate that at high T this energy dependence is a universal function of the the ratio E/T. 
Regarding the level correlations, we trace how the correlator falls off as a function of spatial and 
energy separation between the levels. We also trace how the correlations vanish with increasing 
T. Our most noticeable observation is that for two close energies the correlator changes sign from 
positive (attraction) at small distances to negative (repulsion) at large distances. 

PACS numbers: 73.20.Dx, 71.23.An, 71.55.Jv 
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I. INTRODUCTION 



A distinctive difference between metallic and insulating 
systems lies in statistics of energy levels. It is now a com- 
mon knowledge that in metallic diffusive samples, where 
electrons explore the entire volume in course of their mo- 
tion, energy levels obey Wigner-Dyson statistics (for re- 
view, see Refs. |l]J^). In the deeply insulating regime elec- 
tron states are localized, and can be viewed as classical 
points randomly positioned in space and energy. Thus, 
the levels obey the Poisson statistics. There is presently 
no analytical theory describing how the Wigner-Dyson 
statistics evolve into the Poisson one as the degree of dis- 
order grows. The most interesting critical regime, i.e. the 
level statistics directly at the Anderson transition, was 
first addressed in Refs. H]|, which generated an ongoing 
stream of subsequent studies. In terms of experimentally 
observable quantities, the level statistics correspond to 
the correlation in tunnel conductances at different posi- 
tions of the Fermi level. 

It is also well known that electron-electron interactions 
affect considerably the low-energy behavior of disordered 
systems. In the diffusive regime, they generate a dip in 
the average tunneling density of states at the Fermi sur- 
face. This dip manifests itself as the zero-bias anomaly 
(for review, see Ref. ^|). Interaction effects are even more 
pronounced in the insulating regimea. For Coulomb in- 
teractions between the electrons, the tunneling density of 
states turns to zero as a power law in the vicinity of the 
Fermi level, as argued by Efros and Shklovskifl. This 
phenomenon, dubbed the Coulomb gap, also became a 
subject of a big number of papers, which are too numer- 
ous to be cited here. The works before 1985 are reviewed 
in Refs. §10; recent theoretical activity on the classical 



Coulomb glass (aside from a long-standing problem of 
hopping transport) was_focused on the shape of the den- 
sity of states at zeroE30 and finite temperatups-with a 
finite bandwidth of a bare energy distributionEJO. Re- 
cent studies also addressed the dynamics of formation of 
the Coulomb glasaiiltj. 

Naively, one could think that the effect of interactions 
on the level statistics in the classical case is only to sup- 
press the density of states, which enters as a factor in 
the Poisson distribution. Thus, the resulting distribu- 
tion might be expected to remain Poissonian. In reality, 
however, even in a purely classical system interactions 
cause deviations from the Poisson statistics of energy lev- 
els. To be specific, let us address the following situation. 
Consider a narrow energy strip close to the Fermi level. 
Without interactions, the positions of sites with energies 
within the strip are completely random. The prime ef- 
fect of interactions is an overall depiction of the strip 
(manifestation of the Coulomb gap). Less obvious is the 
question: to what extent, after the energies are modified 
by interactions, the sites within the same strip "feel" each 
other? In other words, how are the fluctuations in their 
density correlated? How does this correlation depend on 
the energy position of the strip? How does it vanish with 
increasing temperature, etc. 

Getting answers to the above questions by means of 
numerical simulations is a much harder task than study- 
ing the shape of the Coulomb gap. In view of this, an 
analytical approach to the problem of level statistics in 
insulators with interactions is vitally important. 

In this paper, we develop an analytical theory of level 
statistics in classical insulators with interactions. 

Such a theory must apparently be based on one of the 
existing analytical treatments of the Coulomb gap prob- 
lem, which we review below. The first approach is due 
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to EfrosQ. The basic underlying assumption of this ap- 
proach is that at zero temperature the ground state of 
the system can be achieved by transpositions of electrons 
within singly occupied pairs of localized states (sites). 
Namely, for an isolated pair of sites with the energies 8\ 
and 82, separated by the distance p\2, the work needed 
to transfer an electron from the occupied site 1 to the 
empty site 2 is 



Ais 



£2 - Si 



(1) 



where V(p) is the pair interaction potential. For the 
ground state, this wotk must be positive for each pair 
of sites (1,2). Efroslij expresses the density of states 
through the probability that all the conditions A12 > 
are fulfilled. The advantage of the approachllj is that it 
yields a closed equation for the density of states. The 
drawback is that it neglects the possibility of reducing 
the total energy by transpositions of electrons within the 
groups of three and more sites.-^An alternative derivation 
of the self-consistent equationta is given in Rcf. |l7|. 

The formalism developed in Ref. [l8] is also restricted 
to transpositions of pairs. The difference between Refs. 
jl6| and [f| lies in the way of pair counting. In Ref. [l^ the 
pairs are ranged according to the distance between the 
constituting sites P12 (referred below as the pair shoul- 
der). The theory describes the evolution of the density of 
states with increasing pi2- Consequently, the final result 
emerges as a solution of a "rate equation" in the limit 
when all the transpositions are carried out (£12 — > 00). 
In contrast to Ref. [16], the formalism of Ref. |l8| preserves 
automatically the total number of occupied and empty 
sites. By this virtue, the latter allows the generalization 
to finite temperatures in a natural way: if the condition 
A12 > is met for a certain pair, the transfer of the 
electron is not forbidden, but happens with a probabil- 
ity exp(— A12/T). Although the formalisms of Refs. [l6| 
and |l8| are different, the results they yield for the zero- 
temperature density of states coincide. 

Another approach was put forward by Johnson and 
KhmelnitskiiE£l. They start from the exact expression for 
the partition function, expand it into the series of Feyn- 
man diagrams, which they subsequently average over dis- 
order term by term. In the one-dimensional case they 
sum up the leading-order logarithmically divergent di- 
agrams. The result of the summation reproduces the 
expressionEHI obtained in the framework of the approach 
of Ref. [l6]. For-two- and three-dimensional cases (D = 
2, 3), the resultsQ for the density of states v(E) oc \E\ D ~ 1 
in the case of the Coulomb interactions (V(p) cx p^ 1 ) 
were reproduced in Ref. 19 by means of an e-expansion 
with e = D — 1. 

Both approaches were designed to describe the average 
density of states. 

Below we study the effect of interactions on the fluc- 
tuations of the density of states in disordered insulators, 
based on the technique developed previously in Ref. [l^. 
In Section [□] we derive the general equations relating the 



average density of states and the level-level correlation 
function for an arbitrary pair interaction potential V(p) 
and finite temperature. Sections [II and IV are devoted 



to the solution of these equations for the case when this 
potential can be treated perturbatively. For this case we 
derive a general formula for the density of states (Sec- 
tion III) with an arbitrary potential and apply it to the 



two particular geometries relevant for experiment: the 
2D disordered insulator located above a planar gate, and 
sandwiched between the two planar gates. The perturba- 
tive expression for level-level correlation function is de- 
rived and analyzed in Section IV. Concluding remarks 
are presented in Section |y|. 



II. GENERAL EQUATIONS 

The main subject of the study of level statistics is the 
level- level correlation function, defined as follows, 



R+(r,E 1 ,E 2 ) = {u + {E x , R)v+(E 2 ,R + r)) c 



(2) 



where the subindex c stands for the irreducible part, 
(AB) C = (AB) - (A)(B), and u+{E,R) is the local den- 
sity of occupied states at the space point R. and energy 
E, 



(3) 



Here we consider a system of electrons which are localized 
at the impurities randomly positioned at space points Ri, 
and having (single-particle) energies rii are the occu- 
pation numbers, assuming the values or 1. It is crucial 
to emphasize that the single-particle energies include 
the interaction-induced shifts from all the surrounding 
sitesO. Note that the definition (Q) differs from the stan- 
dard definition, commonly used in the diffusive regime 
(see e.g. Refs. in two respects. First, only occupied 
states are taken into account. The existence of the Fermi 
level (for which we choose E = 0) breaks the homogene- 
ity in the energy space, and thus Eq. (|2|) depends on 
the two energies E\ and E 2 . Furthermore, the definition 
(^) implicitly assumes that the particles are classical, i.e. 
the spatial extent of the electron wave function is much 
smaller than the average distance between the impurities. 
The level-level correlation function (J2j) depends only on 
the difference r of the space arguments due to the global 
homogeneity. Because of this, the averaging in Eq. (||) 
may be actually viewed as the integration over R. 

To construct the theory, it is convenient to introduce, 
in addition to Eqs. (|^), (0), the density of empty states, 

v- (E, R) = J2^~ ni)S(R - Rt)S (e< - E) , (4) 

i 

and three more correlation functions, 



R-(r,E u E 2 ) = (v-(E u R)v-(E 2 ,R + r)) c 



(5) 
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S(r,E u E 2 ) = (iy+(E 1 ,R)u-(E 2 ,R + r)) c , (6) 



and 



S(r,E 1 ,E 2 ) = {v-{E u R)v + (E 2l R + r)) c , (7) 

They are related by obvious symmetries, 
R-fa-^-Et) = R+(r,E 1 ,E 2 ), and S(r,E u E 2 ) = 
S(r, E 2 , Ei). Thus, the knowledge of the functions R + 
and S is sufficient for the description of the level statis- 
tics. 

Without interactions, the average density of states is 
constant (we denote it with vq), while the fluctuations of 
the local density of states are Poissonian, 

(u ± (E,H)) = Vo f F (±E); (8) 
R ± {r,E u E 2 ) = v Q 8{r)5{E l - E 2 )f F {±E x )f F {±E 2 ); 
S(r,E u E 2 ) = u a 5{r)5{E l - E 2 )f F (E x )[l - f F (E 2 )}. 

Here f F (E) = [exp(E/T) + is the Fermi function. 
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FIG. 1. Transpositions of electrons between single and oc- 
cupied localized states, illustrating the derivation of Eqs. ( p"c| ) 
and (|ll]). The states, characterized by location in energy and 
space, are shown as dashes; occupied states are marked by 
bullets. Processes 1 and 2 are explained in the text. 

In the presence of interactions, Eqs. (@) do not actu- 
ally describe the equilibrium state of the system. To find 
the equilibrium state, we employ the approach of Ref. 
|l8| . Within this approach, the search for the equilibrium 
state is performed by transpositions of electrons within 
pairs of singly-occupied states (i.e. the pairs («,j) with 
rii = 1 and rij = 0). To this end, we introduce an auxil- 
iary variable p (the length of the pair shoulder), and the 
functions g ± (E,R,p). These functions have the mean- 
ing of the local densities of occupied/empty states, which 
are established after the transpositions for all the pairs 
with the shoulders less than p are already carried out. 
For p — the non-interacting values serve as an initial 
condition; the limiting case p — > oo corresponds to the 



equilibrium state of the system. The densities of states 
v (|) and ^ ± are related via 



v ± (E,R) = lim g ± {E,R,p). 

p — >oo 



(9) 



Now we derive the equation for the functions g^. To 
be specific, we consider now the two-dimensional case. 
The derivation is best illustrated by Fig. 1. Consider 
the energy strip (E,E + dE). Imagine that all the 
transpositions with the shoulders shorter than p are al- 
ready performed. As we increase the shoulder from p to 
p + dp, the single-particle energies change due to trans- 
positions, and two processes take place: (i) a certain 
number of occupied sites leave the energy strip (Process 
1 in Fig. 1); (ii) as a result of the transpositions, the 
empty sites with the single-particle energies within the 
strip (E + V(p),E + dE + V(p)) get occupied and arrive 
into the strip (E, E + dE) (Process 2 in Fig. 1). The net 
change of content of the energy strip is responsible for the 
change 5g + (E, R, p) of the density of states. Thus, the 
evolution of g + with p is described by an equation which 
is similar to the rate equation of the kinetic theory, with 
the processes (i) and (ii) corresponding to out-scattering 
and in-scattering terms in the collision integral, respec- 
tively. The actual form of this equation is as follows, 



\ J dn^g+(E,R,p) 
\ dE'g-(E',R + np,p)+ / 

J — oo J E 



dg+(E,R,p) 
dp 

rE+V(p) 



(10) 



dE' 



x g (E ,R + np, p) exp — 



g-(E + V(p),R,p) 



E+V(p) 

E'-E-V(p)\ 



T 



/ dE'g + {E',R + np,p) 
Je 

f E 

/ dE'g + (E',R + np,p)exp 

J OO 



E-E'\ 



T 



where n is a unit vector, and J dn — 2ir. 

Let us establish the correspondence between the r.h.s. 
of Eq. (|l0| ) and the "scattering" processes shown in 
Fig. 1. The first two terms (first square bracket) in the 
r.h.s. describe the scattering out, i.e. departure of levels 
from the energy strip (E, E + dE) due to the transposi- 
tions (Process 1). Indeed, this bracket is multiplied by 
the factor g + , describing the electrons which have been 
transferred. The bracket itself contains the function g~ 
as an integrand, which accounts for empty levels to which 
the electrons are transferred. The argument R + np of 
this function takes into account that the transpositions 
occur within the pairs with the shoulder p. Furthermore, 
in the integration over the energies of all available sites, 
it is taken into account that the transpositions with neg- 
ative work A(E, E') = E' - E - V(p) happen with the 
probability one, while for A(E, E') > the probability 
is given by the Boltzmann factor. The last two terms in 
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the same way describe the scattering in (Process 2). It 
is important that Eq. ( |l0| ) is not averaged over the posi- 
tion R, and thus carries the information about the local 
fluctuations of the density of states. 

Similarly, the "kinetic" equation for g~ takes the form 



dg-(E,R,p) 



dp 



dE'g + {E',R + np,p) 



dn<^g~(E,R : p) 

E-V(p) 



(11) 



E-V(p) 

x g + (E' , R + up, p) exp 
- g+(E-V(p),R,p) 



dE' 



dE'g~ (E',R + np, p) exp - 



E-E' - V(p) 
T 

dE'g-{E',R + np,p) 
E'-E) 



T 



Now we use the equations ( |Io|) and (|Tl]) to obtain 
the coupled equations for the average density of states 
and the level-level correlation functions. The analogy 
between Eqs. (10) and ( pi] ) and the kinetic theory of 
gases suggests the following course of action. We split 
the functions g^ into average (over the space) and fluc- 
tuating parts, 

g ± (E, R, p) = i/± (E) + Sg ± (E, R,p), (12) 

with (Sg^ = 0. First, we average Eq. (|l0|). After this 
procedure, the rate of change of the average value of g + 
is expressed through vj, v~ , and the correlator of the 
fluctuations Sg ± , 

= -up \u+(E) f dE'i>[E + V(p) - E'] V -{E') 



du+(E) 



Op 



- v-[E + V{p)] J dE'ip(E' - £>+(£') (13) 
+ J dE'^[E + V(p)-E']S p (p,E,E') 
dS p (r : E U E 2 ) 



dp 



J dE'ip(E' - E)S p [p, E', E + V(p)} 



where we have introduced the function 
1, 



,B/T 



E>0 

E < 



(14) 



The quantity S p , which enters the r.h.s., is the cross- 
correlator of the fluctuations Sg + and 5g~ , 

S p {r,E 1 ,E 2 ) = (Sg + (E u R,p)Sg-(E 2 ,R + r,p)), (15) 

which, due to the spatial homogeneity, depends only on 
the distance r between the two points. 

Averaging of Eq. ([IT]) yields an analogous equation re- 
lating v~ to vj and the correlator S p . However, we ob- 
serve that the electron-hole symmetry, which applies for 
all average quantities, implies v~(E) = v^{—E). Thus, 
this second equation is unnecessary. 

If we neglect the fluctuation-induced contribution to 
the density of states (terms with S p in the r.h.s.), Eq. 
( |il"| ) turns into a closed equation for iA_ which coincides 
with that obtained previously in Ref. |l§|. 

The equation ( |l3| ) needs to be supplemented by two 
more equations describing the evolution of the correla- 
tion functions S p and R p , defined as 

Rf(r,E 1 ,E 2 ) = (6g ± (E 1 ,R,p)6g ± (E 2 ,R + r,p)). (16) 



To derive these equations, we first subtract from Eq. (|1£ 
for dg + (Ei , R, p) /dp its average, Eq. (13), and call the 
resulting equation A. In the same way, we subtract from 
Eq. ( pd| ) for dg + (E2, R+r, p)/dp its average and call the 
resulting equation B. Following the general approach 
of the kinetic theory of gases, we now multiply A by 
5g~(E 2l R+r 1 p), multiply B by Sg + (Ei,R, p), add them 
up and take the average, disregarding triple correlations 
(see below Section [v]). As a result, we obtain 



-^27rS p (r,E u E 2 ) J dE'^[E x - E' + V{p)]v p (E') - 2*R-[r, E 2 , E 1 + V(p)] j dE'^[E' - E 1 ]u+(E') 

J dn J dE'ip[E\ — E' + V(p)]R p [\r + np\ , E' , E 2 ] 
v-lEk + Vip)] J dn J dE'i/;(E' - Ei)S p [\r + np\,E> \E 2 ] 

2nS p (r, E 1 ,E 2 ) J dE'ip[E' - E 2 + V{p)]v+{E') - 2irR+[r, E 2 - V(p), E 1 ] J dE'^{E 2 - E') Vp {E') 
v p {E 2 ) J dn j dE'^[E' -E 2 + V{p)]R+[\r + np\,E',E 1 ] 

v+[E 2 -V{p)] J dn J dE'^-E^SpWr + np^E^E'A. (17) 
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For the second equation, we ultimately multiply A by and take the average. The resulting equation reads 
dg + (E 2 , R + r, p), multiply B by 8g~{E\,R,p), add up 



dR+(r : E l7 E 2 ) 
dp 

-^2irR+{r, E 1: E 2 ) J dE'i>[E x - E' + V{p)]vj(Ef) - 2nS p [r, E 2 , E x + V{p)\ J dE'^[E' - Ei]u+{E?) 
Vp{Ei) J dn J dE'ijj[Ei — E' + V(p)]S p [\r + np\,E 2 ,E'] 
u-[Ek + V{p)] J dn J dE'i{,(E' - E x )R+[\r + np\,E 2 ,E'} 

2nR+(r, E U E 2 ) J dE '^[E 2 - E 1 + V{p)\v-{E') - 2nS p [r, E U E 2 + V{p)\ J dE'tp(E' - E 2 )v+{E') 
u+(E 2 ) J dn J dE'i>[E 2 -E' + V(p)]S p [\r + np\,E 1 ,E'] 

v-[E 2 + V{p)] J dn j ' dE'^E' - E 2 )R+[\r + np\,E 1) E'}} . (18) 



Equations (|l^), (|i~7|), and ( |l~8|) describe the evolu- 
tion of the functions v p , S p , and i?+ with increasing 
p. Due to the electron-hole symmetry, R~ (r, Ei, E 2 ) = 
Rt(r, —Ei, —E 2 ), and thus the system of the equations 
(|l3|), ([I?]), and ( |l8| ) is closed. For p = 0, when no transpo- 
sitions have yet been performed, the system corresponds 
to the equilibrium state of non-interacting localized elec- 
trons. Thus, Eqs. (||) serve as initial conditions at p = 0, 

v+(E)\ p=o =v f F (E); 



R+(r, E u E 2 )\ p=0 



2nr 



d(r)5(Ei — E 2 )f F (Ei)f F (E 2 ); 



S p (r,E u E 2 )\ p=0 = -^6(r)S(E 1 - E 2 ) 



2nr 

x f F {E x )[l- f F {E 2 )]. 



(19) 



To obtain equilibrium properties of interacting system 
within this approach, one must solve Eqs. (|l3|), (17), 
and (|l8|). The average density of the occupied states 
v + (E) and the correlation functions R + (r, E\ , E 2 ) and 
S(r, Ei, E 2 ) are then obtained as the limit p — > oo of 
the functions u+(E), R+(r,E 1 ,E 2 ), and S p (r,E 1 ,E 2 ), 
respectively. 



III. DENSITY OF STATES 



To this end, we notice that the right-hand sides in Eqs. 
(|j~3|), (p7|), and (18) are effectively small due to the energy 
integration. Thus, in the leading order approximation it 
is sufficient to put the initial conditions ([l9]) into the 
right-hand sides of these equations, and find the devia- 
tions from the initial conditions by integrating them over 



P- 

For the density of states, inserting Eq. ([h]) into the 
right-hand side of Eq. ([l3]) and performing the energy 
integration, we obtain the following equation for , 



Op 



-ivpvlT 

fexp^M / 
[ 1 + expf M 1 + ex P- 
1 



-E-v( P y 



T 



1 + cxp ^ 



exp^P / E" 

l + cx P ^M ln l 1 + eXP T, 

iT^%M ln ( 1 + exp ^ 



(20) 



The cumbersome equations (|l3|), (|l7|), and ([lSl ) contain 
information about the statistical properties of the equi- 
librium state of the system for an arbitrary interaction 
potential V. In this and the next Sections our goal is to 
study weak deviations of the level statistics from the Pois- 
son distribution, and of the average density of states from 
its non-interacting value v$. We solve this problem per- 
turbatively and demonstrate that the desired deviations 
appear already in the first successive approximation. 



valid for an arbitrary temperature0. The total average 
density of states is found as v(E) = v + (E) + v~(E), 
where v + (E) is the solution of Eq. ( |20| ) with the bound- 
ary condition (|l9|), taken for p — > oo. Due to the particle- 
hole symmetry, this becomes v{E) = v + (E) + v + (—E). 

Although Eq. ( ppf ) can not be integrated for an ar- 
bitrary temperature, in two important limiting cases re- 
sults may be obtained in a closed form. 
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A. T = 



2. 2D system with a single gate 



1. General results 



For zero temperature, Eq. (|20|) simplifies to 

-ixpvl [2E + V(p)} 6{-E)6 [E + V(p)] . (21) 



dp 



Integrating Eq. ( j2l| ) , we find for the true average density 
of states 



v(E) = z/ 



rPE 

1 + m> \E\p% -tvuq I dppV(p) 
Jo 



, (22) 



where pe is the solution of the equation V(p) = \E\. 

It is remarkable that although the right-hand side of 
the "kinetic equation" Eq. (|2l]) vanishes after being in- 
tegrated over dE, thus providing the independence of 
the total number of levels on the interactions, the so- 
lution j22] ) does not obey this relation. To realize this, 
we recall that for p — > each realistic potential has a 
Coulomb form, V(p) = e 2 j up. Then it is straightforward 



to check the relation J dE[v(E) — vq] = —7 



1 /2k 2 



in which the right-hand side assumes a finite value due 
to the short-distance behavior of V(p). This seeming 
inconsistency has the following resolution. Let us in- 
troduce a cut-off distance p m i n (minimal separation be- 
tween the sites) which is much smaller than the char- 
acteristic decay length, d, of the potential V(p). Then 
for energies E > e 2 (Kp m i n )~ 1 we have v(E)= v$. With 
finite Pmin the conservation of the number of states, 
J dE[v(E) — vq] = 0, gets restored. This restoration 
occurs due the fact that within a wide interval of high 
energies [e 2 (K/9 m i„) _1 <C E <C e 2 {nd)~ 1 ] the deviation 
v(E) — vq exhibits a shallow maximum. While the height 
of this maximum vanishes with p m i n — > 0, the area under 
the maximum remains finite. 

Conceptually more interesting is the low-energy be- 
havior of the density of states. For E = 0, we obtain 
2/(0) = fo[l — two J °° dppV(p)] < v§. Thus, the aver- 
age density of states at the Fermi level is suppressed by 
interactions. This effect is a precursor of the Coulomb 
gap. 

The perturbative treatment is only valid in the case 
when the correction to the non-perturbed density of 
states z/q, given by Eq. (p2[), is small. Thus, the ex- 
pansion parameter in the theory is vod 2 V(d). If the po- 
tential V(p) falls off at large distances slower than p~ 2 , 
the integral dppV(p) diverges, manifesting that these 
potentials (including the important case of the Coulomb 
potential) are not amenable to the perturbative treat- 
ment at T = 0. 

Now we apply the general result Eq. ( p2| ) to two spe- 
cific realizations of short-range interaction potentials. 



In the presence of a planar gate at a distance d from 
the 2D electron gas, the Coulomb potential is modified 
due to the image charge, 



V(p) 



1 



1 



K \P y/p 2 + Ad 2 



(23) 



where k is the static dielectric constant. In this case the 
decay length is d; for larger distances V(p) falls off as 
p~ 3 . With interaction (|2^) the integration in Eq. ( |2^ ) 
can be carried out in a closed from. However, the ra- 
dius pe is now a solution of the fourth power algebraic 
equation. Then it is convenient to introduce the dimen- 
sionless energy E = 2dnE/e 2 and present the result for 
the density of states as follows, 



v{E) -vq = - 



2nv 2 e 2 d(y(\E\)-l) 2 
k y 2 (\E\) + l 



(24) 



where the function y(E) is the positive solution of the 
equation y 4 + (A/\E\)y 3 -1 = 0. For £ = 0we obtain 
1/(0) — Vq = — 2-Kv^de 2 1 k. In the vicinity of the Fermi 
surface the interaction-induced correction has a singular- 
ity 



v{E) - t/(0) = 2-kvI 



2e 2 d 2 



2/3 



\E\ 1 / 3 . 



(25) 



The energy dependence of the density of states calculated 
from Eq. (p4h is shown in Fig. 2. 




FIG. 2. Density of states ( J24| ) for a 2D system with a single 
gate in units of 2-KV^e 2 dj function of the dimensionless 

energy E at zero temperature. 



3. 2D system between two gates 

We consider now a planar electron system located be- 
tween two infinite planar gates. Denote with d\ and d<i 
the distances to the upper and lower gates. For large 



G 



p 3> (d% + d 2 ) it can be shown that the interaction V(p) 
takes the form 



V(p) 



nK^/2p(d 1 + d 2 ) 

Tip 



: COS 



n(di - d 2 ) 



2{di - 
P 



-(h) 

exp(-/o/a), (26) 



where a = (di + d2)/Tr is the interaction radius; the pa- 
rameter Vq is defined by the second equality in ( pq ) . The 
asymptotic expression (26) is only suitable to describe 
energies \E\ <C Vq. In this case, keeping the two leading 
terms, we write 



Pe = a 



In -r—T In fin — -7 

E 2 V IE 



The leading order in Vp/\E\ terms in the density of states 
are found from Eq. (B2), 



v{E) — Vq = TTVqO 2 



rV2 



-V a 



(27) 



--ln(ln- 



0(1) 



The singularity around the Fermi level is of the type 

v(E)-v(0)ot\E\hi 2 {V /\E\). (28) 

It is remarkable that the singularity survives despite the 
fast decay of the potential V(p) with distance. Moreover, 
since the behavior around E = is determined exactly 
by this long-distance decay, the singularity of the form 
( Pq ) is universal and characteristic of any potential which 
falls off exponentially at large distances. In any realis- 
tic experimental setup the sample is always surrounded 
by a number of gates, and may be characterized by sev- 
eral length scales. The long-p (beyond the longest scale) 
behavior of this potential is always exponential. Thus, 
we conclude that the \E\ In 2 \E\ singularity at small E is 
generic. 



B. High temperatures 

We assume now that the temperature is higher than 
the potential V(p) for any p. In particular, this means 
that the potential is finite everywhere. Expanding Eq. 
( pp| ) in V(p)/T up to the second order (first-order terms 
cancel out) and subsequently integrating over p, we find 
the high-temperature expression for the density of states, 



u(E) = V0 + ^F(E/T), 
where we have introduced the function 



(29) 



F(x) = 



2cosh J (x/2) 



ln(l + e x ) ln(l + e" x ) 



(30) 



and the parameter 
C-- 



d PP V\p). 



(31) 



It is straightforward to check that J dxF{x) — 0, which 
automatically ensures the conservation of the total num- 
ber of states. The effective upper limit in the integral 
( |3l| ) is p ~ d, where d is a characteristic radius of the 
interaction potential. However, the integral diverges log- 
arithmically at small p, since at short distances the in- 
teraction is always Coulomb, V(p) = e 2 j up. In fact, this 



divergence is an artifact of the expansion of (20) with re 



spect to V(p)/T, which is valid only for p^> pr — e 2 /kT. 
Thus, px serves as an ultraviolet cut-off in the integral 
(pl|). Taking this into account, we obtain 



(32) 



6 E/T 




FIG. 3. The universal function F{E/T) (Eq. (|30|)). 

The actual form of the interaction potential affects only 
a numerical factor under the logarithm. So does the con- 
tribution from the short distances p < px- The condi- 
tion of validity of the high-temperature asymptotics (|29| ) 



is T > V{d) - e 2 /nd, 



the argument of logarithm 



m (Jij2j) should be much bigger than one. Another con- 
dition is that the magnitude of the interaction-induced 
correction should be smaller than vq. The latter can be 
expressed as T > v^e^ / n 2 . This condition has a trans- 
parent physical interpretation: T should be bigger than 
the width of the Coulomb gap at zero temperature. 

It is remarkable that the density of states (|2^) depends 
on energy only in combination E/T, and thus the result 
is expressed in terms of the universal function F, shown 
in Fig. 3. The asymptotic expressions for v (E) are 



v(E) = u 



2 r 
ttVqL 

4T 



2 In 2 - 1 



E 2 In 2 
2T 2 



(33) 
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for \E\ < T, and 

v{ E) = V0 + ^e-\ E \l T (34) 

for ^ T. As it is seen from Fig. 3, the maxima of 
v{E) correspond to E = ±2.2T. 

IV. LEVEL STATISTICS 
A. Singular part 

The level-level correlation function is defined by 

K(r,Ei,E 2 ) = (^5(R-R i )6(R + r-R j )8(e i -E 1 )6(E j -E 2 )\ . (36) 

\ ij I c 



K(r,E u E 2 ) = {Av(Ex,K)Ais{E 2 ,TL + r)) 
= lim {S p (r,E 1 ,E 2 ) + S p (r,-E 1 ,-E 2 ) 

p—>oo 

+ R+(r,E 1 ,E 2 ) + R+(r,-E 1 ,-E 2 )}, (35) 

where Av is a fluctuating part of the density of states. 
It is very easy to draw general conclusions about the sin- 
gular part of this expression. Indeed, in terms of the 
localized states wc obtain 



If we now assume that there are no states with different 
energies localized at the same space point (i.e. Ri = Rj 
implies e, = e.j , this assumption effectively serving as an 
ultraviolet cut-off) , it is evident that Eq. ( |36| ) contains a 
singular part, 



K sing (r,E u E 2 ) = y{E 1 )8{E 1 - E 2 ) 



2irr ' 



(37) 



which originates from the term with i = j. On the other 
hand, the remaining terms are regular in r, and can be 
obtained perturbatively from the equations ( |l7j ) and (18), 
to which the initial conditions (119) must be inserted in 
right-hand sides. Below we deal in this way with the 
regular part K reg of the level-level correlation function 



B. Regular part, zero temperature 

For T = 0, the correlation function S reg (r, Ex, E 2 ) de- 
scribes fluctuations of the densities of states at the op- 
posite sides of the Fermi level, E\E 2 < 0; otherwise, it 
equals zero. The result for this quantity is 



S reg {r,E 1 ,E 2 ) = ^9{p* -r). 



(38) 



with the notation p* = P \e 1 \+\e 2 \- In the derivation of 
(pq) we have used the identity 



dn 



S(\r + n P \) 1 
2-k\t + np\ r 



The function R+ eg (r, Ei,E 2 ) is responsible for the corre- 
lations of levels at the same side of the Fermi surface, 



R reg (r,E 1 ,E 2 ) = -fe(p-r), 



(39) 



with p = mhi[pEn Pe 2 ]- Eq. ( |39| ) applies for E\ > 
0, E 2 > as well as E\ < 0, E 2 < 0; therefore we dropped 
the superscript +. 

_The first observation to be made is that the equations 
and (|3^) taken together obey the "detailed balance" 
relation, 



dE x dE 2 K reg (r,E u E 2 ) 



2 dEx dE 2 S reg {r : E u E 2 ) 

JO J-oo 
poo />oo 

2/ dE x I dE 2 R reg (r,E 1 ,E 2 )=0, (40) 
Jo Jo 



for any r. 

Obviously, S reg and R reg are zero without interactions. 
Therefore, they express interaction-induced spatial cor- 
relations between the sites. These correlations are the 
stronger the closer the two energies E\ and E 2 are to the 
Fermi surface: though the amplitudes of the correlators 
are constant, their radii, p or p* , increase as E\ — > and 
E 2 — > 0. The sign of these correlations is negative (level 
repulsion) for the energies of the opposite signs. The 
origin of this repulsion is that the transpositions, which 
build up the ground state of the interacting system, in- 
volve pairs of sites across the Fermi surface. As a result, 
in the ground state, when the transpositions are energet- 
ically unfavorable, the occupied and empty sites tend to 
stay away from each other. 

The positive sign of the correlator R (level attraction) 
can be understood as follows. Close in space occupied 
sites "feel" the repulsion from the same empty sites and, 
hence, tend to gather closer together. Note, that the 
phenomenon of "clustering of like sites" was identified in 
course of numerical simulations with V(p) oc p^ 1 (and, 
thus, fully-developed Coulomb gap) in Ref. |2^. 



8 



C. Regular part, finite temperatures 

The first order corrections to the functions S and R + 
can be actually calculated in a closed form for an arbi- 



trary temperature. Combining them into the level cor- 
relation function according to Eq. (|35|), we obtain an 
expression 



v 2 r 

K reg {r,E 1 ,E 2 ) = -^|-2V>[£i - E 2 + V(r)}f F {E 1 )f F (-E 2 ) - 2ijj[E 2 - E 1 + V{r)]f F {-E 1 )f F (E 2 ) 
+ 1>{Ei - E 2 ) [f F (Ei)f F (-E 2 - V{r)) + f F (-E 2 )f F (E 1 - V(r))} 

+ i>{E 2 - E x ) [f F (E 2 )f F (-E 1 - V{r)) + f F {-E x )f F {E 2 - V(r))}}. (41) 



This is a general result valid for an arbitrary relation 
between the energies E x , E 2 , V(r), and T. It is straight- 
forward to check that for T — Eq. (^Tj) reproduces 
expressions ( |38| ) and (|3^). Also, the detailed balance re- 
lation, 



J dE 1 dE 2 K reg (r,E u E 2 ) = 0, 

is fulfilled for any r and T. 

The expression ( |35| ) for the level-level correlator can be 
presented in a compact form when the energies E\ and 
E 2 coincide, 



K reg {r,E,E) 



v$smh(V(r)/2T 
cosh 2 {E/2T) 



cosh(V(r)/2T\ cosh 2 (E/2T) + sinh(V(r)/2TJ smh 2 {E/2T) 



cosh(£/T) + cosh (V(r)/T 



(42) 



We see that the abrupt cut-off of the correlator (|39|) at 
T = gets smeared at finite temperature. In the vicinity 
of the Fermi level (\E\ <C T) Eq. @ simplifies to 



K reg {r, 0,0) 



■4 , fV(r) 
— tanh ' 



V 2T 



(43) 



The above expression illustrates the decay of the corre- 
lator with increasing either the distance r or the tem- 
perature T. In the opposite limit (\E\ 3> T) the gen- 
eralization of (|3^) to the finite temperatures takes the 
form 



K reg (r,E,E) = 



exp 



(v(r)/T 



- 1 



exp(\E\/T) +2 cosh (V(r)/T 



(44) 



which shows that at finite r, instead of the abrupt termi- 
nation of ( j39l) at \E\ — V(r), the correlator decays with 
departure of E from the Fermi level as exp(—\E\/T). 

Finally, let us consider the behavior of the correlator 
at large distances (V(r) <C T). Expanding ( p2[ ) in the 
first order in V(r)/T yields 



K reg (r,E,E) = 



4Tcosh 2 (S/2T)' 



(45) 



so that [in accordance with Eq. (f43|)] K reg (r, E, E) de- 
cays as V{r). The most important consequence of Eq. 
(E2I) is that at any distances and temperatures the sign 



of the level-level correlator at coinciding energies remains 
positive (level attraction). Below we will see that this is 
not the case for arbitrarily small mismatch E\ — E 2 , To 
illustrate this point we analyze several limiting cases for 
which the general expression (^) for K reg (r, Ex,E 2 ) can 
be simplified. 

Far tail in r, different energies, V(r) <C \E\ — 
E 2 \,T. In this limit we get 



K reg {r,E x ,E 2 ) = - 



16T 2 



exp 



\Ei — E 2 
2T 



1 



cosh(£i/2T) cosii s (E 2 /2T) 
1 

cosh(£ 2 /2T) cosh 3 (^i/2T) 



(46) 



We see that, in contrast to the case of coinciding ener- 
gies, the correlation function is negative (level repulsion) 
and proportional to V 2 (r). 

Far tail in energy, E x = 0, \E 2 \ > T,V(r). In the 
leading order in exp(— \E 2 \/T) we obtain 



K rKq (r,0,E 2 ) 
fV(r) 



v reg \ 
,2 



-yexp 



2T 



2\E 2 \ \ sinh 2 (V(r)/2T) 
cosh(V(r)/2T) ' 



(47) 



This expression indicates that for large enough energy 
separation the repulsion persists for arbitrary relation be- 
tween V(r) and T. 

Crossover from short- to long-distance behav- 
ior. Consider K reg (r, E + A/2, E - A/2) with T, A < 
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\E\. We have seen already that this function (at least for 
A <C T) is positive for low distances r, see Eq. (|43|). We 
are now going to show that, as a function of the separa- 
tion r, this function changes sign at r — pa and becomes 
negative (level repulsion) for long distances. Indeed, for 
Pe -C r < pa (A < y(r) <C |25|) in the leading order we 
find 



K reg (r, E + A/2, 15 - A/2) = i/g exp 
F(r) - A 



A -21251 



2T 



exp 



- 1 



(48) 



which is positive and corresponds to the level attrac- 
tion. The level-level correlation function turns to zero 
for r = pa- For longer distances, the leading order is 
proportional to exp(— 2\E\/T) and becomes (r > pa) 



2-2\E\/T 



K reg (r, E + A/2, 25 - A/2) = 
(l + e~^) (e^)/T _ x ) sinh (m.^ 



(49) 



which matches Eq. (^6|). The levels, indeed, repel each 
other. We emphasize once more that the transition be- 
tween level attraction and level repulsion happens when 
V(r) — A, which corresponds to much longer distances 
than V(r) = \E\. 



D. Level number variance 

We now turn to the question about fluctuations of the 
number of levels within a certain energy strip which was 
formulated in the Introduction. 

For weakly disordered systems, the variance of the 
number of levels in the energy strip of the width E 
serves as quantitative measure of the deviation of the 
level statistics from Poissonian (see e.g. Ref. ^). Pro- 
vided the strip is taken far from the band edges and con- 
tains many levels, this level number variance is a function 
of E only, and does not depend on the absolute position 
of the strip. 

Drawing an analogy with the weakly disordered case, 
we define the quantity E\, E%), which is the variance 
of the number of levels which are located within the spa- 
tial region Q and have energies between E\ and E2 ■ This 
level number variance can be trivially expressed through 
the level-level correlation function, 

E(n,E u E 2 ) 

= [ 2 dEdE' [ dRdR' K(\R — R'\,E, E'). (50) 

J Ei Jfl 

Unlike in weakly disordered case, for many-electron sys- 
tem under study the level number variance depends on 
both Ei and E2 , not merely on the difference E\ — E2 ■ 

Eq. (|5(]) can be significantly simplified if we assume 
that the spatial domain O is sufficiently large, so that for 



relevant energies the correlator K is very small outside 
this domain. Then £ is essentially proportional to the 
volume of f2. Defining the partial level number variance 
(referred below as LNV), a{E u E 2 ) = fi- 1 S(^, Ex, 25 2 ), 
and using the decomposition of K into regular and sin- 
gular parts, we obtain 

a(E u E 2 )~ is {E 2 - -Ei) = <T smg (Ei,E 2 )+a reg (Ei,E 2 ), 

(51) 



where a S i ng and a reg are defined as 



/•E2 

<Jsing{El,E 2 ) = dE\v 
J E x 



(J5i,£ 2 ) 



(E) - vq] , (52) 



dEdE' / 2nrdrK reg {r,E,E'). 



(53) 



I Ei JO 



The value Vo(E 2 — Ei) corresponds to non- interacting 
electrons; a S ing and a reg are the interaction-induced cor- 
rections. The physical meaning of these corrections is 
quite different. Correction a S i ng accounts for the net de- 
pletion of the strip due to reduction of the density of 
states. It does not bear information about deviations 
from the Poissonian statistics. These deviations are de- 
scribed by a reg . 

Generally, both corrections <j S i ng and a reg are compli- 
cated functions of E\, E 2 , T, and the form of the po- 
tential. For this reason, below we consider only several 
instructive particular cases. 

T = 0, symmetric interval. Suppose that the en- 
ergy strip is symmetric with respect to the Fermi level: 
Ei = —25/2, 252 = E/2. At zero temperature the "non- 
Poissonian" contribution to LNV can be evaluated by 
substituting Eqs. Q and © into Eq. M), 



a reg (-E/2,E/2) 



rE/2 r^/'^ r°° 

2nv% / dEi / dE 2 / drr6[p - r] 
Jo Jo Jo 

drr0[p* — r] 



rE/2 



E/2 M poo 

„ , dEi / dE 2 

J-E/2 Jo 



(54) 



With the use of definitions of p(Ei,E 2 ) and p*(Ei E 2 ), 
integration over r and one of the energies in Eq. (p4j) can 
be performed explicitly. Then a reg (—E/2, E/2) takes the 
concise form 



r E 

a reg {~E/2,E/2) = -2-kvI \ dE'[E - E')p%, 

J E/2 



(55) 



This expression is one of the central results of the present 
study. Note that initial formula (|J) for OYe 9 (-£y2, E/2) 



10 



contains two competing terms. The first term is pos- 
itive and originates from attraction between the occu- 
pied as well as between the empty levels. The sec- 
ond term is negative reflecting the repulsion between 
the levels with different occupation numbers. Eq. ( p5| ) 
shows that the second term wins, i.e. the net devia- 
tion from the Poissonian statistics for T = and sym- 
metric energy interval corresponds to the repulsion. To 
proceed further, we choose for concreteness the potential 

/ \ 1/2 

V{p) = Vola/p) exp(— p/a) (p » a), which pertains 

to the planar electron system between two gates. At 
short distances, p -C a, the potential is Coulomb (see 
Section II] ). Consider E <C Vb. Then the integrals ( p2 ) 
and (^||) are determined by distances pe ^> a and can be 
calculated explicitly. The result for LNV reads 



a(-E/2,E/2) - v a E 

- [vo - v(Q)]E + 2tt (In 2-1/2) v%a 2 E 2 In 



(56) 



where the difference i>q — z^(0) is equal to ^l 2 v 2 a 2 Vo/2 as 
follows from Eq. (|2^). It is seen from Eq. ( |56| ) that the 
leading effect of interactions on LNV comes from the de- 
pletion of the strip. The effect of level-level correlations is 
smaller in parameter E/Vq <C 1. Obviously enough, the 
interaction-induced correction to LNV is much smaller 
than i/qE, since vqVqo? «C 1 is the expansion parameter 
of the perturbation theory. 

Note that, in the leading order, a reg (—E/2,E/2) 
is proportional to E 2 ln 2 (Vb/£'). However, this con- 
tribution is canceled by the corresponding part in 
o'sing(— E/2, E/2). As a result, the second term in 
contains ln(Vb/-E) in the first power. 

For E 3> Vo the interaction-induced correction to LNV 
is determined by the short-distance behavior of V(p). It 
depends on ultraviolet cut-off required to satisfy the sum 
rule f[v(E) — i/ ]dE = (see Section III), and is not 
discussed here. 

High temperatures, symmetric interval. The 
above conclusion that correlations cause a negative cor- 
rection to LNV does not hold as the temperature is el- 
evated. To illustrate this, we choose the limit when the 
temperature is much higher than the width of the strip. 
In contrast to the T = case considered above, we 
will also assume that the strip is much wider than the 



crossover energy Vq 



' j na. As it was already men- 



tioned, for T ^> Vo the exact form of the long-distance 
fall-off of the potential appears to be not essential. More 
precisely, the characteristic decay radius a enters the re- 
sults only under the logarithm, as we established in Sec- 
tion |n| while calculating the density of states [see Eq. 
(p2|)]. In this calculation, high temperature introduced 
a short-distance cut-off at p ~ pt = e 2 / kT . The major 
contribution to \v{E) — vq\ came from the wide interval 
of distances pr < p < a within which the potential is 
Coulomb. It appears that the same simplifications are 
valid for the calculation of <j reg {—E/2,E/2). Namely, it 



is safe to use the expansion with respect to V(r)/T in 
the general formula ( fll| ) for the correlator K reg . Sub- 
stituting this expansion into Eq. (|5^) and performing 
integration first over E\ , Ei and then over r we obtain 



,(-E/2,E/2) 



2nv 2 e A E /KaE\ 



(57) 



The positive sign of o~ reg indicates that at high T correla- 
tions cause a super-Poissonian behavior of LNV. The re- 
maining task is to compare Eq. (57) with o~ S ing- The ex- 
pression for <7 S i ng {— E/2, E/2) immediately follows from 
Eq. @. 



,(-E/2,E/2) 



(21n2-l)7reVgE 
4k 2 T 



hi 



KaT\ 



(58) 



Comparing ( p7\ ) and (pq), we see that, besides different 
numerical coefficients, they differ in logarithmic factors, 
which were assumed to be large parameters in the course 
of the calculation. Since T is much bigger than E, the 
logarithmic factor in ( j58| ) dominates. 

In both particular cases considered above, contribution 
of correlations to LNV was parametrically smaller than 
the contribution from the change of the density of states. 
It is interesting to find out whether or not the opposite 
relation is possible within a certain range of parameters. 
Below we illustrate that the opposite relation is indeed 
possible. 

High temperatures, asymmetric interval. Con- 
sider now the case of strongly asymmetric (with respect 
to the Fermi level) energy strip, (E - A/2,E + A/2), 
with E 3> A. The most interesting situation for this 
case occurs when the conditions Vb <C T -C E are 
met. Then [v(E) — v§\ is given by Eq. ( p4| ) and is 
small as cxp(—E/T). Hence, <J s ing is a l so proportional to 
exp(—E/T). As a result, unlike the symmetric-strip case, 
the interaction-induced correction to LNV is dominated 

by <J reg, 

a{E - A/2, E + A/2) - is A 

pE+A/2 poo 

= / dEdE' \ 2irrdrK reg (r,E,E'). (59) 

Je-a/2 Jo 

From Eqs. (g§) and @) we see that the level-level cor- 
relation function for V(r) < E is also exponentially sup- 
pressed. Then the major contribution to ( p9| ) comes from 
the short distances where E < V(r). For such r we can 
replace K reg (r, E±, E 2 ) by K reg (r,E) in Eq. (]5^). More- 
over, for T E < V{p) we can use the zero-temperature 
result ( p9[ ) for level-level correlation function, namely 
K reg (r,E,E) = (vq/2)0(pe — r). Finally, since the cen- 
ter of the strip, E, lies well above Vb, the parameter ps 
is determined by the Coulomb part of the potential, i. e 
Pe = e 2 /kE. Combining all together, we obtain for LNV 



a{E - A/2, E + A/2) = 2/ A 



2E 2 



(60) 



11 



The super-Poissonian result for LNV is not surprising, 
since we have already realized that like sites tend to at- 
tract each other. To estimate the ratio of non-Poissonian 
and Poissonian contributions to LNV, it is convenient to 
rewrite it as (A/ E)(uoe 4 / k 2 E). The first factor is small. 
Less trivial is to realize that the second factor is also 
small. This follows from the condition E ^> T, which 
we exploited, and from the fact that in the perturba- 
tive regime the condition T v$e A / k 2 should hold. As 
it was mentioned above, this condition ensures that the 
Coulomb gap is washed out by temperature. 



V. CONCLUSIONS 

In the present paper we have extended the self- 
consistent description of the disordered interacting sys- 
tem of localized electrons proposed in Ref. [l8] to include 
the correlation of the single-particle levels (both in en- 
ergy and in space). Correspondingly, our results allow to 
incorporate correlations in the calculation of the density 
of low-lying neutral excitations of the system (these exci- 
tations determine the specific heated). We obtained the 
analytical solution of the self-consistent equations for the 
case when, due to the screening by gates, interactions are 
weak enough to be treated perturbatively. Let us discuss 
the conditions when the perturbative description is valid. 

In the absence of screening, the width of the Coulomb 
gap can be estimated by equating the low-energy 
asymptoticsQ v(E) ~ n 2 \E\/e 4 of the density of states to 
the bare density of states vq. This yields \E\ ~ v$e A / n 2 = 
Ec- A perturbative solution applies when the screening 
cuts off the Coulomb interaction at short enough dis- 
tance, d, so that V(d) ~ e 2 /nd 3> Ec- In this case the 
interaction- induced correction to the density of states, 
8v(E), is much smaller than i>q at all energies. In other 
words, the perturbative regime corresponds to a small 
enough uq, which implies a broad distribution of the bare 
energieSr-af sites, i.e. to the limit of a strong external 
disordered. 

The above condition d <C n/i/Qe 2 can be derived from a 
very different consideration. Since the density of states at 
the Fermi level is close to v$ in the perturbative regime, 
the combination Rc — njv^e 2 formally represents the 
static screening radius in two dimensions. Then the con- 
dition of validity of the perturbative description can be 
presented as d <C Rc-, and has a transparent interpre- 
tation that screening is not important for short-range 
interactions. 

At zero temperature and E = we obtained the per- 
turbative result Sv(0)/vq = —ttvq J q dppV(p). Note, 
that the r.h.s. of this expression can be interpreted as 
a probability for an electron at the Fermi level to be 
transposed as the interactions are switched on. Then the 
fact that 5v{Q)/vq <C 1 implies that the portion of sites 
which arc involved in transpositions {i.e. the sites, i, for 
which the transposition with another site, j, is energeti- 



cally favorable) is relatively small. This ensures that the 
deviation from the Poissonian statistics of energy levels 
is weak. 

Obviously, the perturbative description also applies 
when temperature is high enough, T 3> Ec- 

The predictions following from the theory developed, 
which are easiest to test experimentally, concern the en- 
ergy dependence of Sis. The precursor of the Coulomb 
gap should show up in the low-bias behavior of the dif- 
ferential tunneling conductance SG(U) oc 8v(U). For 
tunneling between a metallic gate and two-dimensional 
layer of localized electrons 8G(U) is determined by Eq. 
(§§) (see also Fig. 2) at T = 0, and by F(U/T), where 
the function F is given by ( pCj ) and plotted in Fig. 3, for 
high temperatures. 

Let us discuss the physics left out within the self- 
consistent approach. Firstly, while considering transpo- 
sitions within pairs of sites, we neglected the fact that, 
upon each transposition i — > j, the single-particle ener- 
gies of all other sites change. This change could trigger 
the "secondary" transpositions which were forbidden (at 
T = 0) befop, the transposition i —> j took place (po- 
laronic effected). Secondly, we were concerned only with 
the question whether or not a transposition within a pair 
is energetically favorable. The example given in Ref. |25| 
shows that, even when all pair conditions, A12 > [ see 
Eq. (^) ] are satisfied, the work required for simultane- 
ous transfer of two electrons can be negative (bipolaronic 
effect). Note, however, that both effects are higher or- 
der in vq than the fg-corrections captured by the per- 
turbative solution. For the same reason, we explicitly 
disregarded triple correlations of the type {8gSg8g} when 
deriving Eqs. ( |l7| ) and (|l8|). Terms describing these cor- 
relations are also of the order of and are beyond our 
precision. 

We now give an overview of works discussing spectral 
statistics in interacting systems. One avenue found in 
the literature is to study (mostly numerically) the statis- 
tics of many-particle levels. The early finding was that 
many-particle levels in generic interacting systems obey 
Wigner-Dyson statistics, while integrable ones are Pois- 
sonian distributed!^. Subsequently, statistics of many- 
particle states have been discussed in weakly disordered 
conductorsEH, at the metal-insulator transitionES, and in 
the insulating regimeo, with the general conclusion that 
the many-particle level distribution crosses from the non- 
interacting one to the Poissonian one as the number of in- 
teracting particles increases (thus, the statistics of three- 
particle states are closer to the Poissonian form as those 
of two-particle states). Another opportunity, realized in 
Refs. and is to study the addition spectrum of 
quantum dots in the deeply insulating regime (defined as 
the increment of the many-particle ground state energy 
when electrons are consecutively added to the dot). 

Our results for level-level correlations correspond to yet 
another experimental setup: STM tip moves along the 
surface of the sample and measures the tunneling coa- 
ductance, proportional to the local density of statestH 
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This is unrelated to the addition spectrum of quantum 
dots, since to discuss the addition spectrum, one has to 
account for the shift of the Fermi energy as electrons are 
added. What is more, it is not clear whether the re- 
ally measured addition spectrum of such dots would be 
determined by the ground state energy rather than by en- 
ergies of the single-particle intermediate states, through 
which an electron actually tunnels. What we discuss in 
this paper are also not many-particle states but rather 
quasiparticle excitations (corresponding to the poles of 
the single-particle Green's function); this is an analog 
of the excitation spectrum peaks in the weakly localized 
regimeEH rather than many-body states which compose 
each peak. 

Finally, we remark that recently the issue of the 
Coulomb gap in quantum interactins-systems became a 
subject of several numerical studiesE^TO in connection 
with the quantum Hall effect and metal-insulator tran- 
sition in two dimensions. The results of our paper are 
strictly classical and are valid in the deeply insulating 
regime, i.e we assumed that the overlap integral between 
the localized wave functions is smaller than the smallest 
energy scale in the problem; however, to match them to 
the zero-bias anomaly in the weak disorder regime oac 
inevitably has to pass the quantum insulating domainEi 
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